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ABSTRACT 

SN 2002es is a peculiar subluminous Type la supernova (SN la) with a combination of observed 
characteristics never before seen in a SN la. At maximum light, SN 2002es shares spectroscopic prop- 
erties with the underluminous SN 1991bg subclass of SNe la, but with substantially lower expansion 
velocities (~6000 km s~ x ) more typical of the peculiar SN 2002cx subclass. One month after maxi- 
mum light, spectra of SN 2002es show low velocities as seen in SN 2002cx. Photometrically, SN 2002es 
differs from both SN 1991bg-like and SN 2002cx-like supernovae. Although at maximum light it is sub- 
luminous (M B = -17.78 mag), SN 2002es has a relatively broad light curve (Am 15 (B) = 1.28 ± 0.04 
mag), making it a significant outlier in the light-curve width vs. luminosity relationship (the Phillips 
relation). From Arnett's law, we estimate a 56 Ni mass of 0.17±0.05 M Q synthesized in the explosion, 
relatively low for a SN la. One month after maximum light, we find an unexpected plummet in the 
bolomctric luminosity. The late-time decay of the light curves is inconsistent with our estimated 56 Ni 
mass, indicating that either the light curve was not completely powered by 56 Ni decay or the ejecta 
became optically thin to 7-rays within a month after maximum light. The host galaxy is classified 
as an SO galaxy with little to no star formation, indicating the progenitor of SN 2002es is likely from 
an old stellar population. We also present a less extensive dataset for SN 1999bh, an object which 
shares similar photometric and spectroscopic properties. Both objects were found as part of the Lick 
Observatory Supernova Search, allowing us to estimate that these objects should account for 2.5% of 
SNe la within a fixed volume. We investigate theoretical models to probe the nature of SN 2002es, 
but find that current models are unable to explain all of its characteristics. 
Subject headings: supernovae: general - supernovae: individual (SN 2002es, SN 1999bh) 



1. INTRODUCTION 

Type la supernovae (SNe la) are the runaway ther- 
monuclear explosions of carbon-oxygen white dwarfs. 
SNe la are characterized spectroscopically by an absence 
of hydrogen and the presence of intermediate-mass el- 
ements (e.g., silicon, sulfur, oxygen) and iro n-group el- 
ements (iron, cobalt) fe.g.. iFilippenkol L1997I . and refer- 
ences therein) . A majority of spectroscopically identified 
SNe la form a class of objects with a standardizable lu- 
minosity, allowing for their use as distance indicators. 
Application of SNe la on extragalactic scales led to the 
disco yery that the Universe is accelerating in i ts expan- 
sion (jRiess et all 119981: IPerlmutter et all I1999D . Subse- 
quent app lication of large samples of SNe la out to high 
redshifts dWood-Vasev et all 120071: iHicken et al.l P2 009a 



Kessler et al. 2009: Amanullah et al.ll2010: Sullivan et alj 



2011at ISuzuki et alj 12012ft has led to precise estimates 
of cosmological parameters when combined with mea- 
surements of baryon acoustic oscillations (BAO) and 
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anisotropy in the cosmic microwave background (CMB). 

The cosmological application of SNe la is predicated 
on the relationship between the peak absolute magni- 
tude of a SN, t he width of its light curve, and its color. 
iPhillipsI (| 1993( 1 found that SNe with slowly declining 
light curves had a larger luminosity at maximum 
light. Applying corrections for light-curve width and 
SN color (as a measurement of host-galaxy extinction 
and intrinsic scatter in SN colors) has allowed SNe la 
to be ac curate distance indicators to within 10% in 
dista nce (jJha et al J 120071: IGuv et all 120071: iConlev et alj 
2008). Ther e are indica ti ons t h at including s pectral 



information (iFoley et al.1 f2 008: 
Bailev et al.l l2009t iFolev fc Kasenl 
201 It iSilverman et al. 



Wang et all 



20111 iBlondin et al.1 



2009a; 



formation (iKellv et al.l 



2012a) and host-gala xy 

mm 



ax y m- 
ISullivan et all l20ToT: 
iLampeitl et al.1 l2010f) further improves distance esti- 
mates. 

Despite the ability to standardize the luminosity of 
SNe la based on observed light-curve properties, a sig- 
nificant fraction of SN la events have peculiar char- 
acteristics, including s ome overluminou s and underlu- 
mi nous objects (e.g ., [Filippcn kol 119971). In p articu- 
lar . iFilippenko et all (|i992D andlLcibundg ut et all (|1993Q 
found that optical light curves of SN 1991bg evolved 
rapidly, and its peak luminosity was ~2 mag fainter than 
that of normal objects. Moreover, the maximum-light 
spectrum of SN 1991bg showed strong Ti II absorption, 
indicating a relatively cool photosphere, and the expan- 
sion velocity at maximum light as measured from the 
absorption minimum in the blueshifted Si II A6355 fea- 
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ture was ^10,000 km s , slightly lower than expan- 
sion velocities measured for normal SNe la (^11,000 - 
13,000 km s" 1 ) (FilippcnkoeLaU [HIS) . Since the ini- 
tial identification of SN 1991bg as a subclass of SNe la, 
many members belonging to the subclass have been iden- 
tified by the SN community. Studies of the host-galaxy 



-type galaxies ( 


Howellll2001t 


ILi et all 12011a 


), leading to 



iGane shaling am et al.l H( 
the suggestion that the progenitors of SN 1991bg-like 
SNe come from old stellar populations. There is currently 
debate about the cosmological utility of SN 1991bg- 
like objects as standardizable candles (| Jha et al.l 120071 : 
IGuv et al]l2007h . 

More recently, a range of propertie s for peculiar sublu- 
minous SNe la have been discovered. iFolev et al.l ()2010bf ) 
presented evidence that SN 2006bt spectroscopically re- 
sembled SN 1991bg, but photometrically resembled a 
normal SN la. SN 2006bt was discovered at a pro- 
jected distance of 33.7 kpc from the nucleus of its early- 
type host galaxy . The Palomar Transient Factory (PTF; 
lLaw et alj|2009t ) has published data on two peculiar sub- 
luminous objects, both found at large distances from the 
likely host galaxy. PTF 09dav was an abnormally sublu- 
minous SN la (Mb = —15.44 mag) with a narrow light 
curve fo und 41 kpc from its host galaxy (jSullivan et al.l 
l2011bH . iMaguire et"aT1 (|2011| ) presented data on the sub- 
luminous PTF lOops which shared many similarities with 
SN 2006bt. PTF lOops had a broad light curve and was 
found at a projected distance of 148 kpc from the nom- 
inal host. All objects had spectral features that match 
those of SN 1991bg, although PTF 09dav had particu- 
larly slow expansion velocities of 6100 km s" 1 . 

Objects sim ilar to SN 2002cx (jFilippenkol 120031: 
ILi et al.l l2003bf) form another subclass of subluminous 
peculiar SNe la. These objects have maximum-light 
spectra similar to those of overluminous objects like 
SN 1991T, characterized by weak Si II A6355 features 
and dominated by Fe III lines indicating a hot photo- 
sphere. However, the expansion velocities of these ob- 
jects at maximum light are ~6000 km s , indicating an 
explosion with low kinetic energy per unit mass. There 
appears to be a great diversity among SN 2002cx-like 
objects, with a distribution of absolute luminosity and 
kinet ic energy (jNaravan et al.1 120111 : iMcClelland et all 
2010). SN 2008ha is the faintest member of the subclass, 
with M v = -14.2 mag a nd velocities of ~ 4000-5000 
km s _1 at maximum light (|Folev et al.ll2010"al ). 

Here we report our observations of SN 2002es, an ob- 
ject somewhat spectroscopically similar to SN 1991bg, 
further adding to the puz zle of sublum inous SNe la. 
SN 2002es was discovered (|Li et al.l 120021) on unfiltered 
CCD images at ~16.3 mag as part of the Lick Observa- 
tory Supernova Search (LOSS) with the 0.76- m Katzman 
Automatic Imaging Telescope (KAIT; iFilippenko et al.1 
[2001 : ILi et aTll2003aft on 2002 Aug. 23.5 (UT dates are 
used throughout this paper) in UGC 2708. Its J2000 co- 
ordinates are a = 03 h 23 m 47 s .23 and <5 = +40°33 m 53.5 s , 
which is 19" W and 26" N of the galaxy nucleus 
(|Li et al.l I2002D . Subsequent o ptical spectrosco p ic ob - 
serv ations on 2002 S e p. 3 by iChornock et al.l (2002) 
and iMatheson et all (12002ft classified SN 2 002es as a 
SN 1991bg-like object (jFilippenko et al.lll992t ) near max- 
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Fig. 1.— KAIT V-band image of SN 2002es. The field is 6.7' X 
6.7'. The SN and comparison stars are marked. The labels of the 
comparison stars correspond to the numbers in Table [T] 

imum light based on the presence of strong Ti II, Si II, 
and O II absorption features, but with an expansion ve- 
locity of ~6000 km s _1 as measured from the absorp- 
tion minimum of the Si II A6355 feature. This is notably 
lower than typical expansion velocities of ^11,000 km s _1 
found in normal SNe la and SN 1991bg -like objects an d 
more typical of SN 2002cx-like objects (|Li et al.ll2003rJ ). 

Noting the peculiarity of SN 2002es, our group started 
a photometric and spectroscopic campaign to document 
the evolution of this unique object. We present BVRI 
photometry for SN 2002es starting a week before max- 
imum light collect ed as part of the LOSS SN la pho- 
tometry program (jGaneshalingam eTall [2010h . Late- 
time photometry was obtained with the 2.3-m Bok tele- 
scope at Steward Observatory on Kitt Peak in Arizona 
to constrain the late-time decay of the light curve. We 
also present an extensive optical spectral series covering 
the evolution of SN 2002es from maximum light to two 
months after maximum. In addition, we present evidence 
that SN 1999bh is a SN 2002es-like event. 

The rest of the paper is structured as follows. Our data 
and reduction techniques are presented in fj^l with a de- 
tailed analysis of our photometry and spectroscopy in £j3j 
In Sj4l we discuss possible physical interpretations for the 
observed properties of SN 2002es and compare SN 2002es 
to theoretical models. We summarize our results in f}5] 

2. OBSERVATIONS 

2.1. Photometry 

Broadband photometric follow-up observations of 
SN 2002es started on 2002 Aug. 24.48 (one week before 
B max ) with KAIT in the BVRI bands. SN 2002es was 
monitored with a 1-2 day cadence for the first month af- 
ter discovery, resulting in well-sampled light curves. Ad- 
ditional BVRI photometry from the 1-m Nickel Tele- 
scope at Lick Observatory was obtained to complement 
the KAIT data. Late-time data in BVR were taken us- 
ing the 2.3-m Bok telescope. A l/-band image of the field 
from KAIT is shown in Figure [1] 
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Images were bias subtracted and flat fielded using stan- 
dard procedures in IRAF0 SN 2002es was sufficiently 
close to its host galaxy that galaxy subtraction was re- 
quired to disentangle SN light from galaxy light. Tem- 
plates for the field were obtained for the KAIT, Nickel, 
and Bok telescopes after SN 2002es faded beyond the 
detection limit of the telescope (>2 yr after discovery) 
and subtracted from the data images. The flux from the 
SN was measured in comparison to local field stars us- 
ing point-spread function (PSF) fitting photometry. The 
uncertainties in our photometric measurements were es- 
timated by randomly injecting artificial stars with the 
same PSF and magnitude as the SN into the data im- 
ages. The images with artificial stars were rerun through 
our galaxy-subtraction and photometry routines, and the 
scatter in doing this 20 times was adopted as the photo- 
metric measurement error. 

Instrumental ma gnitudes wer e color -corrected to the 
Landolt system (|Landoltl 119921 |2009[ ) using the aver- 
age color terms measured from multiple photometric 
nights. The magnitudes of local field stars were cal- 
ibrated against Landolt standards on six photometric 
nights. The results are reported in Table [T] The final 
photometric error is the photometric measurement error 
added in quadrature with the error in the calibration of 
our field stars. The final photometry is presented in Ta- 
bleU 

In addition to the photometry presented here, we sup- 
plement our light curves with UB VRI data presented for 
SN 20 02es in the CfA3 photometry sample (iHicken et al.l 
2009b) taken with the 1.2-m telescope at Fred Lawrence 
Whipple Observatory (FLWO) in Arizona operated by 
the Harvard Smithsonian Center for Astrophysics (CfA). 
The agreement between the KAIT and CfA3 data is 
within ~ 0.05 mag in all bands, with no evidence of a 
systematic offset. The light curves from all telescopes 
are presented in Figure [5J 

2.2. Spectroscopy 

Low-resolution spectra of SN 2002es were obtained 
using the Kast dual spectrograph mou nted on the 3- 
m Sh ane Telescope at Lick observatory ([Miller fe Stond 
1993), the FAST spectrograph mounted on the 1.5-m 
Tillinghast telescope at FLWO dFabricant et all ]1998[) 



and LRIS on the 10-m Keck I telescope ( Oke et all 



1995). 



All of our observations were taken at the opt imal paral- 
lactic angle to minimize differential light loss ( Filippcnko 
fl98l . 

All spectra of SN 2002es were re duced using stan- 
dard CCD processing techniques (e.g.. lFolev et al.ll2003l ; 
IMatheson et al.ll2008f) . Processing and extraction of the 
one-dimensional spectrum were perf ormed i n I RAF using 
the optimal extraction algorithm of IHorn c (1986). The 
wavelength calibration was obtained by fitting low-order 
polynomials to calibration-lamp spectra. Our spectra 
were flux calibrated using our own IDL routines. Cor- 
rections for telluric absorption fe atures were made us- 
ing spectrophotometric standards (jWade &i Hornd[T988f) 

9 IRAF: The Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astron- 
omy (AURA), Inc., under cooperative agreement with the National 
Science Foundation (NSF). 
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Fig. 2.— UBVRI light curves of SN 2002es. Data from the 
KAIT, Nickel, and B ok telescopes are plotted as circles. Data from 
Hicken ct al. (2009b) arc plotted as diamonds. For comparison, we 
plot the SN 2002cx-like SN 2005hk (solid line), the SN 1991bg-like 
SN 1999by (dashed), and the normal SN 2005cf (dot-dashed line). 
Comparison light curves have been shifted to have the same peak 
magnitude and phase as SN 2002es. Note the particularly fast 
decline in all bands at t > +30 d compared to the other objects. 

taken at roughly the same airmass as the SN observa- 
tion. Table |3] presents a summary of our spectroscopic 
observations. 

2.3. Host Galaxy 

The NASA/IPAC Extragalactic Database (NED)E3 
lists two discordant values for the heliocentric redshift 
of UGC 2708 from different sources. The NED web- 
page for UGC 2708 lists z hc iio = 0.028 as the de- 
fault redshift determin ed from marginal measurement s 
of the H I 21-cm line bv lMonnier Ragaig ne et al.l (2003). 
NED warns under the "E ssential Note" section that 
Ide Vaucouleurs ef a l. (1991) find Zhciio = 0.018 from op- 
tical lines. We obtained a spectrum of UGC 2708 with 
Kast at Lick Observatory on 2006 July 21.4 to determine 
the actual redshift of the host galaxy. We find a helio- 
centric redshift of Zhci = 0.0182 ± 0.0001 from measure- 
ments of weak, narrow Ha + [N II] emission lines, i n goo d 
agreement with the value of lde Vaucouleurs et alj {1991). 
The heliocentric redshift of UGC 2708 corresponds to 
zcmb = 0.0177 (in the frame of the cosmic microwave 
background). We adopt an uncertainty of 300 km s _1 
to account for any peculiar motions induced by gravi- 
tational interactions with neighboring galaxies. For a 
standard ACDM cosmology with Q m = 0.27, f2 A = 0-73, 
w = —1, and Ho = 73.8 km s _1 Mpc -1 (iRiess et alJ 

10 http://ned.ipac.caltech.edu/ . 
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TABLE 1 

Photometry of Local Standard Stars 



ID 


a (J2000) 


<5 (J2000) 


B(a B ) (mag) 


V(a v ) (mag) 


R{<tr) (mag) 


I (a!) (mag) 


-^calib 


1 


03:23:57.81 


+40:36:08.0 


17.403 


(Oil) 


16.497 


(008) 


15.940 


(012) 


15.462 


(007) 


4 


2 


03:23:35.94 


+40:36:00.3 


16.470 


(010) 


15.675 


(007) 


15.206 


(007) 


14.744 


(012) 


5 


3 


03:23:59.75 


+40:36:02.8 


17.817 


(010) 


17.174 


(Oil) 


16.785 


(012) 


16.331 


(014) 


3 


1 


03:23:52.38 


+40:35:25.2 


17.440 


(010) 


16.844 


(007) 


16.482 


(Oil) 


16.104 


(Oil) 


1 


5 


03:23:39.24 


+40:35:16.1 


14.293 


(009) 


12.933 


(014) 


12.191 


(002) 


11.545 


(Oil) 


2 


6 


03:23:42.58 


+40:35:14.0 


16.442 


(007) 


15.702 


(003) 


15.265 


(005) 


14.851 


(Oil) 


5 


7 


03:23:41.02 


+40:35:12.2 


17.978 


(011) 


17.314 


(Oil) 


16.905 


(013) 


16.482 


(010) 


5 


8 


03:23:50.00 


+40:35:08.6 


15.777 


(011) 


15.119 


(005) 


14.727 


(007) 


14.325 


(Oil) 


4 


9 


03:23:40.51 


+40:34:56.4 


18.618 


(013) 


17.970 


(004) 


17.579 


(010) 


17.035 


(008) 


2 


10 


03:23:51.59 


+40:34:50.8 


18.657 


(013) 


17.786 


(006) 


17.292 


(008) 


16.865 


(006) 


3 


11 


03:23:54.37 


+40:34:26.8 


14.523 


(012) 


14.040 


(Oil) 


13.736 


(013) 


13.407 


(006) 


4 


12 


03:23:58.96 


+40:34:29.2 


17.496 


(012) 


16.721 


(Oil) 


16.254 


(013) 


15.735 


(013) 


4 


13 


03:23:51.57 


+40:34:03.6 


16.938 


(010) 


16.180 


(005) 


15.744 


(Oil) 


15.296 


(015) 


4 


11 


03:23:37.25 


+40:33:52.3 


17.115 


(009) 


16.126 


(006) 


15.542 


(007) 


15.047 


(012) 


4 


15 


03:23:56.96 


+40:33:44.5 


15.783 


(012) 


15.258 


(012) 


14.905 


(007) 


14.554 


(010) 


4 


1(5 


03:23:58.62 


+40:33:41.4 


16.993 


(010) 


16.280 


(Oil) 


15.833 


(008) 


15.395 


(015) 


4 


17 


03:23:39.26 


+40:33:28.8 


16.795 


(010) 


16.222 


(005) 


15.863 


(008) 


15.503 


(012) 


4 


18 


03:23:58.71 


+40:33:28.4 


17.935 


(Oil) 


17.021 


(003) 


16.477 


(008) 


16.003 


(007) 


4 


19 


03:24:01.27 


4-40-33-24 7 


18.038 


(013) 


16.979 


(010) 


16.402 


(008) 


15.814 


(017) 


3 


20 


03:23:37.30 


+40:33:17.8 


16.652 


(013) 


15.842 


(Oil) 


15.360 


(009) 


14.935 


(013) 


4 


21 


03:23:56.13 


+40:33:03.5 


18.860 


(010) 


17.906 


(005) 


17.374 


(012) 


16.956 


(008) 


3 


22 


03:23:42.20 


+40:32:48.4 


15.821 


(Oil) 


14.705 


(003) 


14.090 


(005) 


13.512 


(009) 


5 


23 


03:24:00.85 


+40:32:47.5 


17.769 


(009) 


17.064 


(005) 


16.658 


(006) 


16.233 


(005) 


2 


24 


03:23:56.43 


+40:32:43.2 


17.724 


(Oil) 


16.895 


(004) 


16.377 


(009) 


15.836 


(013) 


5 


25 


03:23:49.70 


+40:32:38.5 


16.919 


(009) 


15.795 


(004) 


15.101 


(007) 


14.540 


(003) 


i 


26 


03:23:51.07 


+40:31:60.0 


16.680 


(Oil) 


15.987 


(005) 


15.583 


(007) 


15.170 


(013) 


5 


27 


03:23:41.64 


+40:31:46.7 


18.034 


(012) 


16.877 


(Oil) 


16.157 


(012) 


15.481 


(014) 


2 


28 


03:23:37.03 


+40:31:41.7 


14.570 


(007) 


13.860 


(007) 


13.427 


(012) 


13.003 


(016) 


3 


29 


03:23:55.50 


+40:31:41.6 


18.669 


(Oil) 


17.703 


(012) 


17.062 


(Oil) 


16.456 


(007) 


4 


30 


03:23:38.90 


+40:30:56.3 


17.954 


(013) 


17.261 


(010) 


16.836 


(013) 


16.417 


(014) 


4 



Note. — la uncertainties (in units of 0.001 mag) are listed in parentheses. 



|2011| ), we find a luminosity distance — 73.20 ± 4.41 
Mpc and a distance modulus fi = 34.32 ± 0.12 mag. 
SN 2002es exploded at a projected distance of 11 kpc 
from the nucleus of UGC 2708. 

The original IAUCs th at spectroscopically c lassi 
SN 2002es (jChornock et al.1 [20021 : iMatheson et all[200" 
predate the source for the incorrect redshift listed on 
NED and likely used the c orrect redshift supplied by 
Ide Vaucouleurs et al.l ([1991). The correct redshift for 
SN 2002es was also used in the cosmology analysis of 
iHicken et ail (|2009aft . 

UGC 2708 was included in a s tudy of SN la host-galaxy 
properties by iNeill et al . (2009). The authors estimated 
galaxy properties by fitting template-galaxy spectral en- 
ergy distributions (SEDs) to multi-wavelength photom- 
etry from the Sloan Digital Sky Survey (SDSS) and the 
Galaxy Evolution Explorer (GALEX). For UGC 2708, 
the authors estimated a minimal amount of active star 
formation and E(B — U)host = 0.0 mag. This finding, 
along with the absence of Na I D absorption at the red- 
shift of the host galaxy in spectra of SN 2002es, indicate 
that the reddening due to the host galaxy is negligible. 

UGC 2708 was also observed spectroscopically as part 
of SDSS. Using the publicly available line-fiux measure- 
ments, we find an [N II] /Ha ratio consistent with that 
of a low-ionizat ion nuclear emission-line region (LINER; 
iHeckmanl I1980D or a composite galaxy. Coupled with 
the results from INeill etldl (|2009l) . UGC 2708 is likely a 
LINER with no active star formation. 

3. RESULTS 
3.1. Photometry 



The light curves of SN 2002es are displayed in Figure 
El along with light curves of the normal S N la 2005cf 
fWang et al.l l2008t iGaneshalingam et al.l I2010D the 
SN 2002cx-like SN la 2005hk (iPhillips et al.TT2007h. and 
the S N 1991bg-like SN la 1999bv (jGaneshalingam et al.1 
[2010l) . 

Basic photometric properties for the light curves of 
SN 2002es are reported in Table |U All values were mea- 
sured using a fifth-order polynomial fit directly to the 
data. Uncertainties were estimated using a Monte Carlo 
routine to produce 50 realizations of our dataset. Each 
individual dataset realization was produced by randomly 
perturbing each photometry data point using its pho- 
tometric uncertainty assuming a Gaussian distribution. 
Light-curve properties were measured for each realiza- 
tion, and the final measurements were found by taking 
the mean and standard deviation of the set of simulated 
realizations. 

The light curves of SN 2002es share some character- 
istics in common with the subluminous SN 1991bg sub- 
type. SN 2002 es lacks a prominent shoulder in the R an d 
/ bands (jFilippenko et al.lll992tlLeibundgut et alj fl993). 
The secondary maximum often found in R and / (and 
more prominently in the near-infrared) is attributed to 
the cooling of the ejecta to temperatures where the tran- 
sition from Fe III to Fe II becomes favorable, redistribut- 
ing flux from shorter wavelengths to longer wavelengths 
(jKasenl [20061. This transition occurs earlier in cooler 
SNe. From m odels of the radiative transfer within SNe, 
IKasenl (|2006f) finds that the timing and strength of the 
shoulder is dependent on the distribution and amount 
of 56 Ni within the ejecta. Models with a completely ho- 
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TABLE 2 
Photometry of SN 2002es 



JD - 2,452,000 B (mag) V (mag) R (mag) / (mag) Telescope 



510 


98 


17 


682 


(033) 


17 


245 


(028) 


17.007 


(025) 


16 


835 


(031) 


KAIT 


511 


98 


17 


625 


(029) 


17 


148 


(027) 


16.882 


(022) 


16 


711 


(029) 


KAIT 


512 


97 


17 


498 


(033) 


17 


050 


(028) 


16.786 


(027) 


16 


622 


(030) 


KAIT 


514 


00 


17 


482 


(034) 


17 


009 


(029) 


16.711 


(029) 


16 


521 


(033) 


KAIT 


515 


01 


17 


453 


(035) 


16 


918 


(026) 


16.669 


(020) 


16 


470 


(030) 


KAIT 


516 


01 








16 


752 


(111) 


16.552 


(097) 


16 


375 


(033) 


KAIT 


517 


02 


17 


372 


(104) 






16.587 


(055) 


16 


347 


(046) 


KAIT 


518 


01 


17 


307 


(044) 


16 


778 


(045) 


16.465 


(090) 


16 


374 


(022) 


KAIT 


519 


02 


17 


373 


(020) 


16 


817 


(024) 


16.503 


(020) 


16 


332 


(026) 


KAIT 


520 


01 


17 


369 


(020) 


16 


762 


(020) 


16.493 


(022) 


16 


318 


(026) 


KAIT 


520 


99 


17 


414 


(026) 


16 


790 


(020) 


16.453 


(020) 


16 


303 


(026) 


KAIT 


521 


99 








16 


829 


(031) 






16 


211 


(111) 


KAIT 


523 


02 


17 


506 


(023) 


16 


829 


(021) 


16.484 


(020) 


16 


320 


(031) 


KAIT 


524 


02 


17 


550 


(026) 


16 


844 


(020) 


16.498 


(020) 


16 


339 


(037) 


KAIT 


524 


95 


17 


729 


(020) 


16 


925 


(020) 


16.514 


(020) 


16 


311 


(020) 


1 m 


525 


95 


17 


805 


(020) 


16 


985 


(020) 


16.539 


(020) 


16 


336 


(020) 


1 m 


528 


97 


18 


094 


(020) 


17 


165 


(020) 


16.687 


(020) 


16 


400 


(020) 


1 m 


529 


00 


18 


142 


(030) 


17 


164 


(028) 


16.640 


(023) 


16 


400 


(023) 


KAIT 


530 


01 


18 


255 


(040) 


17 


193 


(025) 


16.746 


(020) 


16 


446 


(023) 


KAIT 


531 


01 


18 


323 


(039) 


17 


293 


(031) 


16.774 


(020) 


16 


492 


(022) 


KAIT 


531 


99 


18 


500 


(049) 


17 


354 


(020) 


16.814 


(020) 


16 


533 


(029) 


KAIT 


533 


uo 


18 


583 


(071) 


1 / 


525 


(080) 


16.914 


(105) 


16 


588 


[101) 




535 


00 


18 


779 


(037) 


17 


556 


(031) 


16.977 


(033) 


16 


621 


(025) 


KAIT 


536 


02 


18 


824 


(053) 


17 


629 


(032) 


17.042 


(033) 


16 


641 


(028) 


KAIT 


537 


99 


18 


876 


(094) 


17 


717 


(050) 


17.122 


(030) 


16 


732 


(043) 


KAIT 


540 


87 


19 


041 


(155) 


17 


836 


(063) 


17.248 


(043) 


16 


792 


(058) 


KAIT 


547 


98 


19 


429 


(167) 


18 


275 


(084) 


17.715 


(036) 


17 


110 


(055) 


KAIT 


550 


97 


19 


610 


(082) 


18 


464 


(061) 


17.910 


(032) 


17 


346 


(033) 


KAIT 


551 


99 


19 


744 


(032) 


18 


570 


(025) 


17.964 


(020) 


17 


417 


(042) 


1 m 


552 


84 


19 


771 


(043) 


18 


631 


(027) 


18.010 


(021) 


17 


505 


(025) 


1 m 


553 


96 


19 


834 


(100) 


18 


602 


(070) 


18.116 


(038) 


17 


615 


(058) 


KAIT 


561 


91 


















18 


481 


(334) 


KAIT 


562 


91 








19 


447 


(118) 


19.167 


(109) 


18 


508 


(123) 


KAIT 


566 


87 








19 


579 


(161) 


19.445 


(125) 


18 


961 


(120) 


KAIT 


573 


86 








20 


155 


(171) 












KAIT 


577 


92 








20 


811 


(307) 












KAIT 


592 


89 


23 


180 (150) 


22 


711 


(110) 


22.34 (110) 








Bok 



Note. — 1<t uncertainties (in units of 0.001 mag) are listed in parentheses. 



TABLE 3 

Log of Optical Spectral Observations for SN 2002es 



UT Date 


Phase a (d) 


Telescope/Instrument 


Exposure Time (s) 


Observer 15 


2002 Sep. 03.5 


+3 


Lick/Kast 


600 


AF, RC, BS 


2002 Sep. 05.5 


+5 


FLWO/FAST 


1200 


MC 


2002 Sep. 06.5 


+6 


FLWO/FAST 


1200 


MC 


2002 Sep. 10.5 


+10 


FLWO/FAST 


1200 


PB 


2002 Sep. 12.5 


+12 


FLWO/FAST 


1200 


MC 


2002 Sep. 13.3 


+13 


Lick/Kast 


1200 


RF, SJ, MP 


2002 Sep. 28.4 


+28 


FLWO/FAST 


1200 


MC 


2002 Oct. 01.3 


+30 


Lick/Kast 


1800 


AF,RF 


2002 Oct. 08 


+37 


Keck I/LRIS 


400 


AF,RC 


2002 Nov. 08 


+67 


Keck I/LRIS 


1800 


AF, RC,SJ, BB 


2002 Nov. 11 


+70 


Keck I/LRIS 


1800 


AF, RC 



a Rest-frame days relative to the date of B max , 2002 Aug. 31.8 (JD 2,452,518.3), rounded to 
the nearest day. 

b AF = A. Filippenko, BB = B. Barris, BS = B. Swift, MC = M. Calkins, MP = M. Papcnkova, 
PB = P. Berlind, RC = R. Chornock, RF = R. Foley, SJ = S. Jha 
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mogenized composition and with a small amount of 56 Ni 
result in an /-band light curve with no discernible sec- 
ondary peak or shoulder. Instead, the two peaks merge 
to produce a single broad peak. Given the similar abso- 
lute magnitudes of SN 2002es and SN 1991bg, the lack of 
a shoulder or secondary maximum in the R and I bands 
may be a consequence of low 56 Ni production, if the light 
curve is powered by the decay of 56 Ni. 

The timing of maximum light in each band is similar to 
that of SN 1991bg-like SNe. In normal SNe, maximum 
light in the I band precedes that of B by a few days 
(as evident in the light cur ves of SN 2005cf in Figure [5]). 



iTaub cnbcrgcr et al 



()2008l ) find that for the SN 1991bg- 
like SN 2005bl, peak brightness in UBVRI occurred in 
successive order of bluest to reddest filter with the date of 
maximum light in each band separated by ~1 d, similar 
to what is seen in SN 2002es. 

The B-band peak magnitude of SN 2002es is 17.33 ± 
0.02 mag after correcting for a color excess of E(B — V) = 
0.183 mag from Milky Way extinction (Schlegel et al. 
1998). We assume negligible host-galaxy extinction 
based on the classification of UGC 2708 as an SO galaxy 
and no evidence of Na I D absorption at the host-galaxy 
rcdshift. This corresponds to an absolute magnitude of 
Mb = —17.78 ± 0.12 mag, comparable to other sublumi- 
nous SNe la, but brighter than SN 1991bg-like objects 
(jTaubenberger et al.ll2008h . The absolute magnitude of 
SN 2002es is brighter than that of the SN 1991bg-like 
SN 2005bl by -0.5-1 mag in BVRI . 

Despite these similarities with SN 1991bg-like SNe la, 
the light curves of SN 2002es are significantly broader 
in all bands compared to those of most subluminous ob- 
jects. Using the decline in magnitudes between maxi- 
mum light and 15 days after maximum light in the B 
band as a proxy for light-curve width ([Phillips! Il993f ). 
we measure Ami^(B) — 1.28 ± 0.04 mag. Other 
SN 1991bg-like SNe typica lly have Ami 5 (5) « 1.9 mag 
(jTaubenberger et al.ll2008h . 

The light curves of SN 2002es also share similarities to 
the SN 2002cx-like SN 2005hk. Both objects are sublumi- 
nous events compared to normal SNe la with light curves 
that are broader than those of SN 1991bg-like objects and 
lack a shoulder in the R and / bands. SN 2002es has a 
slower B-band decline than most SN 2002cx-like objects. 
SN 2005hk had Am 15 (B) = 1.56 ± 0.09 mag. 

As the outer ejecta turn transparent, the light curve is 
expected to be powered by the thermalization of 7-rays 
produced by the decay of 56 Co (—50 d after explosion). 
This is typically ob served as a linear fading in all bands. 
Lcibundgut (2000) find typical decay rates for SNe la 
of 0.014 mag d" 1 in B, 0.028 mag d^ 1 in V, and 0.042 
mag d _1 in /. For SN 1991bg, the B-band decline is 
marginally faster at 0.019 mag d _1 and slower in / at 
0.040 mag d^ 1 . The U-band decline has been found to be 
fairly constant am ong normal and SN 1991bg-like SNe la 
(|Leibundgutl[2000f ). 

Using data at t > +30 d (relative to maximum light in 
B band), we measure a decline of 0.040 ± 0.004 mag d _1 
in B, 0.081 ± 0.004 mag d" 1 in V, 0.101 ± 0.004 mag d" 1 
in R, and 0.099 ± 0.004 mag d _1 in /. These rates are 
substantially faster than expected for an object powered 
by the decay of 56 Co, even after accounting for the de- 
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Fig. 3. — M B as a fun ction of Ami 5 (B) for 76 SNe from 
Ganeshalingam et al. (2010) corrected for host-galxay extinction 
using values derived from the MLCS distance-fitting software. 
Overplotted a s a so lid line is the quadratic Phillips relation from 
Phillips ct al. (1999) shifted along the ordinate to match the data. 
The dashed lines correspond to the lcr scatter about the relation- 
ship. SN 2002es and SN 1999bh (red triangles) are clear outliers 
in the relationship. For comparison, we include other subluminous 
peculiar objects SN 2005hk, PTF 09dav, and PTF lOops. 

homologously, and cast doubt on whether SN 2 002es is 
necessarily a thermonuclear explosion. In §4.3[ we dis- 
cuss possible explanations for such a fast decline. 

In Figured we plot Am 15 (B) vs. M B for SN 2002es 
along with 76 SNe taken from iGaneshalingam et al.l 
(2010) using host-galaxy extinction values deter- 
mined using the Multicolor Ligh t-Curve Shape method 
(MLCS2k2.v006;[jh^iian[2003). The light-curve widt h 
vs. luminosity relationship from [Phillips et all (1999), 
adjusted to Ho = 73.8 km s^ 1 Mpc^ 1 , is overplotted as a 
solid line with la scatter about the relation indicated by 
dashed lines. We also i nclude measuremen ts for the pe- 
culiar SNe SN 2005hk (IPhillips et a l.l l200"l. PTF 09dav 
(iSullivan et al.l l2011b[ ). and PTF lOops (jMaguire et al.l 
120111 ). SN2"002es is an obvious outlier (5a) in the Phillips 
relation which predicts that SN 2002es should be —1.3 
mag brighter than the observed peak magnitude. The 
amount of host-galaxy extinction required to explain this 
discrepancy is unlikely given the absence of Na I D ab- 
sorption in SN 2002es spectra at the redshift of the host 
galaxy. 

3.2. Color Curves 

In Figure [4] we plot the color evolution of SN 2002es 
along with SN 2005cf, SN 2005hk, and SN 1999by for 
comparison. All objects have been corre cted for Milky 
Way extinction using the dust maps of Schlegel ct al. 
(1998). We adopt a host-galaxy co lor excess of E(B — 
V) host = 0.10 mag for SN 2005cf (I Wang et al.ll2009bD. 
■E(£ -U)wt = 0.10 mag for SN 2005hk (IChornock et all 
I2001I. and E(B - V) host = 0.0 mag for SN 1999by 
( Garnavich eTall I2001 . We assume no host-galaxy ex- 
tinction for SN 2002es. 

The B — V color curve of SN 2002es does not match 
that of any of our comparison objects. The color evo- 
lution is most similar to that of the SN 1991bg-like 
SN 1999by, but with significant differences in B — V af- 
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TABLE 4 

Photometric Properties of SN 2002es 



Filter 


JD of max - 2,452,000 


Mag at max a 


Peak abs. mag b 


Amis (mag) 


B 


518.30 ±0.28 


17.33 ±0.02 


-17.78 ±0.12 


1.28 ±0.04 


V 


519.21 ± 0.28 


16.78 ±0.02 


-18.15 ±0.12 


0.74 ± 0.02 


R 


520.96 ± 0.30 


16.47 ± 0.02 


-18.35 ±0.12 


0.57 ±0.02 


I 


521.39 ±0.33 


16.29 ± 0.02 


-18.40 ±0.12 


0.37 ±0.02 



a Not corrected for Milky Way or host-galaxy extinction. 

b Corrected for Milky Way extinction and assuming no host-galaxy extinction. 
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Fig. 4. — Optical color curves of SN 2002es corrected for E(B — 
V) = 0.183 mag from Milky Way extinction. Solid circles are data 
take n with the KAIT and Nickel telescopes and solid diamonds are 
from lHicken et a l. (2009b). Shown for comparison are SN 2005hk, 
SN 1999by, and SN 2005cf corrected for extinction using the red- 
dening values provided in their respective reference. SN 2002es is 
significantly redder in all colors compared to SN 2005cf. None 
of the comparison color curves provides an adequate match to 
SN 2002es. 

ter maximum light in B. Before t(B max ), SN 2002es 
and SN 1999by share a similar B — V color evolution 
that is considerably redder than that of SN 2005cf and 
SN 2005hk. However, at ^5 d past maximum light in 
B, SN 1999by quickly becomes redder in B — V while 
the color evolution of SN 2002es is much more gradual. 
The V — R and V —I color evolution for SN 2002es and 
SN 1999by appear much more similar. There is evidence 
that SN 2002es becomes bluer at t > +35 d, although 
the data are noisy. 

3.3. Spectral Properties 

In Figure [SJ we show our time series of spectra covering 
the evolution of SN 2002es starting 3 d after maximum 
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Fig. 5. — Montage of spectra of SN 2002es. All spectra have 
been shifted to the rest frame by the recession velocity of the host 
galaxy and corrected for E(B — V)mw = 0.183 mag. The phase 
relative to maximum light in the B band is indicated to the left of 
each spectrum. 

light in the B band and extending out to two months af- 
ter maximum light. The absence of Ha and the presence 
of Si II A6355 in early-time spe ctra identify SN 2002es as 
a SN la (e.g.. lFilippenkdll997i ). Strong Ti II absorption 
around 4200 A is commonly associated with the sublumi- 
nous SN 1991bg-like subclass of SNe la (Filippcnk o et all 

In Figure [6J we show our earliest spectrum 
of SN 2002es compared to spe ctra of the nor- 
mal SN 200 5cf flWang et all [20 09b). the subluminous 
SN 1991bg (Turatto e t all 119961 via the online SUS- 
PECT databastQ) and the SN 2002cx-like SN 2005hk 
(|Phillips et al.1 120071 ). At this phase, SN 2002es shares 
the most similarities with SN 1991bg. However, the ex- 

11 http://suspect.nhn.ou.edu/~suspect . 



Ganeshalingam et al. 



1 1 1 

- ivy\ 

Ti 

- *ywv 


S C/Co? 


■ i 1 

^05hk (+4) 
? Ca 


- 

05cf(-l) - 










: y 

— . i 


1 


1 


^9ib g (+2) : 


4000 


6000 


8000 


10000 



Rest Wavelength (A) 

Fig. 6. — Spectrum of SN 2002es taken +3 d after maximum 
light in B, compare d to the spectroscopically norma l SN la 2005cf 
llWang et al.ll2009bl l. the SN 2002cx-like S N 2005hk (IPhillips et al.l 
1200711 . andthe subluminous SN 1991bg HTuratto et al.l 11999 ) at 
similar phases. All objects have been shifted to the rest frame by 
their host-galaxy recession velocity and corrected for both Milky 
Way and host-galaxy extinction using the values in their respec- 
tive references. Major spectroscopic features are identified. The 
spectrum of SN 2002es shares many similarities with SN 1991bg, 
except the former has significantly slower ejecta velocities. In par- 
ticular, both objects show strong O I and Ti II features. The ejecta 
velocities of SN 2002es are more similar to those of SN 2005hk. 




4000 6000 8000 o 10000 

Rest Wavelength (A) 

Fig. 7. — Spectrum of SN 2002es taken + 13 d after maximum 
light. For c omparison we show SN 2005cf (Wang et al. 2009b), 
SN 2 005hk ( IPhillips et al.ll2007f) , and SN 1991bg ITuratto et al.l 
1996) at a similar phase indicated in parentheses to the right of 
each spectrum. All spectra have been corrected for both Milky 
Way and host-galaxy extinction using the values in their respective 
references. Major spectroscopic features are identified. 

pansion velocity of the photosphere as measured from 
the minimum in the blueshift of the P-Cyngni profile 
of Si II A6355 is 6000 km s _1 , significantly lower than 
typical values of ~ 11,000 km s" 1 fo r normal SNe la 
(|Wang et alJ I2009at iFolev et all l20TTI) and SN 1991bg - 
like SNe of -10,000 km s" 1 (jTaubenberger et all [20081 ) . 
SN 2002 cx-like objects have low ejecta veloc ities of —6000 
kms^ 1 (|Li et al.ll2003bl : IPhillips et al.ll2007[ ). comparable 
to what is found in SN 2002es. 




4000 6000 8000 o 10000 

Rest Wavelength (A) 

Fig. 8.— Spectrum of SN 2002es taken +37 d after max- 
imum light using LRIS on Keck. For c omparison we show 
SN 2005cf (Wan g et al.M2009bl). SN 2005hk UPhillips et al.|[2007l ). 
and SN 1991bg l|Turatto et**aTTl996l) at a similar phase indicated in 
parentheses. All spectra have been corrected for both Milky Way 
and host-galaxy extinction using the values in their respective ref- 
erences. Major spectroscopic features are identified. 

The spectrum of SN 2005hk only covers 3600-7400 A. 
Within this spectral range, it shows rather similar fea- 
tures to those of SN 2002es, but with notable differences 
as well. In particular, it has a weaker Si II A6355 ab- 
sorption feature, and no obvious Ti II absorption trough 
near 4200 A. 

About two weeks after maximum light, SN 2002es be- 
gins to show more similarities with SN 2002cx-like ob- 
jects. In Figure [71 we show our spectrum of SN 2002es 
from 13 d after maximum light compared to those of 
other objects at a similar phase. At this phase, Si II 
A6355 absorption is harder to discern. The Si II/ Na I 
complex around 5700 A is weaker than what is seen in 
SN 1991bg. Similar to SN 2002cx-like objects, Fe II 
AA4555, 5129 absorption lines become more apparent and 
the Ca II H&K lines are similar in strength. However, 
O I A7774 and the Ca II near-infrared (NIR) triplet are 
stronger in SN 2002es compared to SN 2005hk. 

In Figure HI we show our +37 d Keck spectrum along 
with comparison objects. We begin to see emission dom- 
inate in the case of the NIR Ca II triplet. The low ex- 
pansion velocities in SN 2002es and SN 2005hk make it 
possible to see narrow features which are usually hard to 
discern due to the smearing of broad high-velocity fea- 
tures in normal SNe la. The spectrum is dominated by 
permitted iron-group elements, but line blending makes 
it difficult to uniquely identify features. 

We show our combined +67 d and +70 d spectrum of 
SN 2002es in Figure |H1 along with spectra of the nor- 
mal SN 1994D, SN 2005hk, and SN 1991bg at compa- 
rable phases. SN 2002es continues to share the most 
similarities with SN 1991bg. Despite the rapid evolu- 
tion of the light curves, the spectra do not show signs 
of being completely nebular. We continue to see con- 
tinuum emission and absorption in our combined spec- 
trum. There is a hint of forbidden emission lines, pos- 
sibly [Ca II] or [Fe II] in the region around 7200 A, but 
we do not detect other prominent forbidden iron emis- 
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Fig. 9.— Combined spectrum of SN 2002es from +67 d and +70 
d after maximum light. For comparison we show t he normal SN 
1994D (from our datab ase of spectra), SN 2005hk ijPhillips et al.1 
I2007T) , and SN 1991bg IITuratto et al.| [l996l at a similar phase in- 
dicated in parentheses to the right of each spectrum. All spectra 
have been corrected for both Milky Way and host-galaxy extinc- 
tion using the values in their respective references. At this phase, 
SN 2002es looks similar to SN 1991bg, but with narrower features. 
Major spectroscopic features are identified. 
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Fig. 10.— Spectra of SN 2002es at +3 d, + 13 d, and +37 d 
after maximum light compared to our best SYNOW fit for each 
epoch in red. Ions that contribute to major features are labeled 
according to the blueshifted absorption minimum of that feature. 
Our tentative identification of C II is indicated by a question mark. 



sion features common l y seen in nebular SN la spectra 
(jStanishev et all l2007t iLeloudas et~aTl 12009ft . We note 
that the permitted Ca II emission is not particularly- 
strong, unlike PTF 09dav and the class of Ca-rich ob- 



jects (Kasli wal et al.H201lD . 



3.4. SYNOW Modeling 

The supernova spectral synthesis code SYNOW 
(jFisher et al.l ll~997) is useful in identifying the different 
ion contributions in SN spectra. We use SYNOW as a 
tool to analyze our +3 d, +13 d, and +37 d spectra to 
determine the elemental composition of the photosphere. 
Our synthetic spectra compared to our observed spectra, 
along with identification of major spectral features, can 
be found in Figure [TU1 

For our +3 d spectrum, we set the blackbody temper- 
ature to 13,000 K with a photosphcric velocity of 6200 
km s . Many of the ions present in our SYNOW spec- 
trum are commonly found in the spectra of SN 1991bg- 
like objects: Ca II, T i II, O I, Si II, S II, and Na I. 
Sulliv an et al.l (|201 lbf ) find evidence of Sc II in spectra 
of PTF 09dav, a peculiar subluminous SN la similar in 
some characteristics to SN 2002es. We find that includ- 
ing Sc II does not improve the fit to the spectrum of 
SN 2002es. 

Just redward of Si II A6355 is a notch that looks 
suspiciously like C II A6580. The presence of carbon 
in pre-maximum spectra has recently be en found in a 
substantial fraction (~20-30%) of SNe la dParrent et al.l 
[20111: iThonias et alj|201lt iFolatelli et al.ll2012ft at veloc 
ities slightly higher than those of intermediate-mass el- 
ements such as Si and Ca. If the notch is due solely to 
C II, it has a velocity of 4200 km s~ l , much lower than 
the velocities of other ions. With SYNOW, we can fit this 
feature reasonably well with a blend of C II and Co II 
with a velocity of 6200 km s — 1 . We do not see evidence of 
C II A7234. The notch is also present in our spectrum at 
+13 d, but it is not generally seen in normal SNe la that 
exhibit carbon at early times before maximum light. Per- 
sistent carbon features lasting two w eeks were seen in the 
excep tionally luminous SN 2009dc ()Taubenberger et all 
120111 ). However, given the subluminous nature of 
SN 2002es, it shares few characteristics in common with 
SN 2009dc (lYamanaka et al.ll2009t ISilverman et a,l.H2011t 
iTaubenberger et al.l 1201 1[) and similar overly luminous 
SNe la (IHowelll 120011: IHicken et al.l [20071 IScalzo et al.l 
[20irl . 

In an a nalysis of the pre -maximum spectrum of 
SN 2006bt, iFolev et al.l (|2010bft found possible evidence 
for C II at slower velocities (^5200 km s _1 ) than other 
ions (^12,500 km s _1 ). However, unlike the case in 
SN 2002es, they found that the feature disappeared after 
maximum light. This feature, whether attributed to C II 
or not, hints at a possible connection between the two 
transients. 

In our +13 d spectrum, we set the blackbody tem- 
perature to 7500 K with a photospheric velocity of 5500 
km s . The low velocities of the ions allow for the de- 
tection of a forest of absorption lines normally smeared 
out by higher velocities. We see less evidence for the 
presence of Si II and S II. Absorption from iron-group el- 
ements such as Fe II and Ti II becomes more prominent. 

By +37 d, we use a blackbody temperature of 5900 
K with a photospheric velocity of 3500 km s . At this 
phase, we begin to see the SN become nebular and signif- 
icant line blending at short wavelengths from iron-group 
elements. We still detect continuum emission, indicating 
that the SN is not completely nebular, but we begin to 
see emission from calcium. 



10 



Ganeshalingam et al. 



0.8 




0.8 1.0 1.2 1.4 1.6 1.8 2.0 
A m 15 (B) [mag] 



Fig. 11. — Top panel: The silicon ratio near maximum light, 
1Z(Si II), versus Ami5(B). Blue diamonds are low- velocity objects 
, red squares are high- vel ocity objects, and y ellow triangles are 
faint objects as defined by Benctti ct al. (2005). SN 2002es joins 
PTF lOops and SN 2006bt as outliers in the relationship between 
7£(Si II) and Ami 5 (B). SN 2002es has a 7£(Si II) more typical of 
a faint object, despite having a broad light curve. Bottom panel: 
Gradient in the velocity of the Si II feature versus Ami5(B). Based 
on the clustering of this plot, SN 2002es would be classified as a 
HVG object despite being subluminous. 

3.5. Quantitative Measurements 

SNe la can be broadly classified by a cl ustering anal- 
ysis of spectral and phot ometric features ([Benetti et al.l 
120051: Branch et al.l [20061) . Most quantitative measure- 
ments make use of S i II fe atures at 5972 A and Si II 
A6355. INueent et all (|1995[ ) introduced ft(Si), the ratio 
of the depth of Si II A5972 to A6355 near maximum light, 
which is found to correlate with light-curve width. Using 
lZ (Si) and the ve l ocity gradient of the Si II A6355 feature, 
v, IBenetti et all (120051 ) found that SNe la could be bro- 
ken into three distinct subclasses: low velocity gradient 
(LVG), high velocity gradient (HVG), and FAINT. 

We estimate 7^(Si) = 0.55 ± 0.05 from our earli- 
est spectrum of SN 2002es taken +3 d after maximum 
light. In the top panel of Figure ITT1 we plot Arn^jB) 
against 7^(Si) for the objects in IBenetti eTaTl (2005) 
along with SN 2002es and other peculiar subluminous 
objects. While SN 2002es has a 7^(Si) value similar to 
that of FAINT objects, its measured Ami^(B) does not 
match those of FAINT objects. Compared to other pe- 
culiar sublumi nous SNe la, SN 200 2es is most similar 
to PTF lOops (jMaeuire et all 1201 ID . both of which are 
significant outliers in the relationship between Ami^B) 
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Fig. 12. — Bolometric luminosity of SN 2002es (dashed line) 
constructed using BVRI data and bolometric corrections from 
SN 1991bg. The shaded region indicates the ltr errors. Plotted 
for compariso n is the bolometric luminosity of SN 2005cf (longer 
dashed line; I Wang eTaTl l2009bT) and SN 1999by (dot-dash line; 
Garnavich ct al. 2004). The solid blue and dotted r ed c urves are 
models for the radioactive luminosity using Equation 1101 assuming 
different 7-ray trapping efficiencies with an initial Mji; = 0.17 Mq. 
The red dotted curve assumes full trapping of 7-rays producing a 
late-time slope that follows the cobalt decay rate. The blue curve 
is our best fit. Neither curve can reproduce the steep late-time 
decay in the bolometric light curve of SN 2002es. 

and n(Si) found by IBenetti et~aTI pOOl . 

We measure v = 98 ± 5 km s _1 d _1 (statistical error 
only) from 6 spectra taken in the range +3 < t < +13 d. 
In the bottom panel of Figure QTJ we plot Am^(B) 
agains t v for SN 2002es and the objects in IBenetti et al.l 
(2005), coded by subclass. Based on our measured v 
and Am 15 (B), SN 2002es falls into the HVG group. 
FAINT objects also have high velocity gradients, but 
they also have a larger Ami5(B) (i.e., a narrow light 
curve). SN 2002es would fall into the FAINT classifica- 
tion if it had a narrower light curve. Both panels indicate 
that SN 2002es has a unique combination of photometric 
and spectroscopi c properties that lie outside the classifi- 
cation scheme of IBenetti et al.l (|2005f) . 

4. DISCUSSION 
4.1. Bolometric Luminosity 

Estimating the bolometric luminosity is particularly 
difficult without observations in the ultraviolet (UV) and 
IR, which contribute a sizable fraction of energy for nor- 
mal SNe la. Unfor tunately, [/-band photometry from 
iHicken et al.l ()2009bf ) starts ~5 days after i?-band max- 
imum light. We only make use of BVRI data t o con- 
struct a "quasi-bolometric" (jNomoto et al.|[T990[ ) lumi- 
nosity light curve for SN 2002es. 

We estimate the quasi-bolometric luminosity by warp- 
ing spectra representing the SED of SN 2002es corrected 
for Milky Way extinction to match our mu lti-color pho- 
tomet ry using a method similar to that of iHowell et al.l 
(2009). We assume there is no host-galaxy extinction 
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(see H2.3p . For each epoch of BVRI photometry, we 
warp the spectrum using a third-order spline with knots 
at the effective wavelength of each filter to match the 
photometric colors. In instances where we are missing 
photometry from one band, we do a linear interpolation 
between the nearest photometric epochs. We note that 
there is a fairly large gap in our i?-band photometry 
at +40 < t < +70 d, making our interpolation uncer- 
tain. The resulting warped spectrum is integrated over 
the range 4000-8800 A (i.e., from the blue limit of the B 
band to the red limit of the / band) to obtain the optical 
flux for the photometric epoch. The flux is converted to 
a luminosity using the distance reported in H2.3I 

Ideally, a spectral series of SN 2002es matched with 
each photometric epoch (as opposed to a single spec- 
trum) should be used for the most accurate results. The 
spectra from FAST presented in this paper do not ex- 
tend to the I band, and thus are not useful for estimat- 
ing the bolometric luminosity. Given the spectroscopic 
similarities between SN 1991bg and SN 2002es, we used 
th e above proce d ure w ith the SN 1991bg spectral series 
of iNueent et all (|2002f ) artificially redshifted by 0.01 to 
match the slower expansion velocities of SN 2002es. 

A bolometric correction is required to turn our quasi- 
bolometric luminosity based on optical data into a bolo- 
metric luminosity that accounts for energy emitted in 
the UV and IR. Our SN 1991bg spectral series covers 
1000-25,000 A and provides a reasonable first-order ap- 
proximation of the SED of SN 2002es. We calculate the 
fraction of flux emitted in the range 4000-8800 A com- 
pared to the total integrated flux for each spectrum in 
the series. This provides us with an estimate of the flux 
missed by only using optical data to construct our quasi- 
bolometric light curve as a function of phase. We ap- 
ply this correction to arrive at our final bolometric light 
curve. At maximum light, the BVRI data account for 
— 70% of the total flux. Given the uncertainties in our 
estimates, we include a 20% systematic error in our error 
budget for the bolometric luminosity. 

We estimate L bol = (4.0 ± 0.9) x 10 42 erg s' 1 at -1 d 
after maximum light in the B band. This is slightly 
larger than t he bolometric luminosity of other sublumi- 
nous SNe la (jTaubenberger et al.ll2008ft and -25% of the 
luminosity of SN 2005cf (|Wang et al.H2009bD . 

The luminosity of thermonuclear SNe is powered by 
the energy deposition of 7-rays and positions produced 
by the radioactive decay chain 56 Ni -> 56 Co -> 56 Fe. At 
maximum light, the rate of energy deposition into the 
expanding ejecta is roughly eq uivalent to th e luminosity 
of the SN (i.e., Arnett's l aw; lArnett] 11982ft . Following 
IStritzinger fc Leibundgutl (]2005l ). we can write Arnett's 
law as 

L hol = a x (6.45 e -'^( 8 - 8d ) + 1.45 e - t ^ m ' 3d )) x (1) 

jxlO ergs , 

where a is a correction factor of order unity to Arnett's 
law and t r is the time between explosion and maximum 
light (i.e., the bolometric rise time). 

Assuming the luminosity of SN 2002es is powered by 
the decay of 56 Ni, we can estimate the amount of 56 Ni 
synthesized in the explosion. We set a = 1. Unfortu- 



nately, we do not have tight constraints on the date of 
explosion. Prior to the first detection of SN 2002es on 
2002 Aug. 23.2 (—9 d before the time of maximum L^ \), 
KAIT obtained an unfiltered image of the field on 2002 
Aug. 12.5 with a limiting magnit ude of 18.5 (|Li et al.1 
120021 ) . iGaneshalingam et all (|2011h found that the sub- 
luminous SN 1999by had a rise time to maximum light in 
B of 13.33+0.40 d, consistent with ot her estimates for the 
rise tim e of SN 1991bg-like SNe from lTaubenberger et al.l 
(2008), while a typical, normal SN la has a rise time of 
~18 d. The S-band light curve of SN 2002es has a slower 
rise to maximum li ght than SN 1999by a nd matches the 
rise of SN 2005hk. lPhiiliDS et al.l ((20071) measure a rise 
to £?-band max imum of 15 ± 1 d fo r the SN 2002cx- 
like SN 2005hk. iMaguire etall (|201lD constrain the rise 
time of PTF lOops, an object that shares similarities to 
SN 2002es, to ~19 d. Based on the unique spectroscopic 
and photometric peculiarities of SN 2002es, it is unclear 
which object serves as the best guide to determining the 
rise time. As a compromise between the different possible 
rise times, we adopt a rise time to bolometric maximum 
of t r = 16 ± 3 d. 

From Equation [U we estimate -Mm = 0.17 ± 0.05 M Q 
synthesized in the explosion which falls at the low 
en d of the range 0. 05 < M Ni < 0.87 M found 
by IStritzinger et al.l (|2006bl ) for a sample of SNe la. 
SN 1991bg and SN 1999by synthesized -0.1 M Q of 56 Ni. 

4.2. Energetics 

We can estimate the ejected mass, M e j, and the 
kinetic energy, Eq, of the explos ion usi n g the rise 
time. Following the treatments of lArnettl (|1982f) and 
iPinto fc Eastaianl (|2000al lbl). we have t 2 . oc nM e -Jv 8 , 
where k is the mean opacity and v s is the ejecta ve- 
locity. If we compare SN 2002es to a normal SN la with 
M e j w 1.4 Mg, o, = 10 4 km s" 1 , and t d = 18 d, and 
assume they have similar opacities, we have 

M ^°- 66MQ (lTd) 2 ( 6OOOkms-0 - ® 

The largest uncertainty is our calculation is from our 
estimate of the rise time of t r = 16 ± 3 d. Including this 
uncertainty, we find M c j = 0.66 ± 0.25 for our nominal 
values, significantly lower than the canonical values for 
a SN la. 

We can then calculate the kinetic energy 
£k .o 4xl05 . erg (^) 3 (_^) 3 . <„ 

Including the uncertainty in our rise time gives E^ = 
(2.4 ± 0.9) x 10 50 erg for our nominal values. 

The estimated ejected mass and kinetic energy are sig- 
nificantly lower than the canonical values for a SN la in 
part due to the slow expansion velocities measured from 
the Si feature. We again caution that our estimate for 
the rise time is not well constrained and can range from 
13 d to 19 d. 

Foll owing similar arguments m ade by IHowell et al.l 
42I||) and ISilverman etaU (|20TTI) . we can place con- 
straints on the white dwarf (WD) progenitor mass, Mwd, 
based on the energetics of the explosion. The amount of 
energy produced by burning carbon and oxygen up to 



12 



Ganeshalingam et al. 



intermediate-mass elements (IMEs) and iron-group ele- 
ments (IGEs) must equal the energy required to unbind 
the WD and the kinetic energy of the ejecta: E nuc = 
Ek + -^binding- The energy released from fusing a WD 
consisting of equal parts carbon and oxygen to iron-group 
elements is 1.55 x 10 51 erg M ^ 1 , and fusion up to 28 Si 
releases 76% as much energy (|Branchlll992f) . Assuming 
the entire WD is burned in the thermonuclear explosion, 
the fraction of iron-group elements, /ige, and the frac- 
tion of intermediate- mass elements, /ime, will add to 1. 
The energy released from nuclear fusion can be written 
as 



E nuc = 1.55 x 10 51 erg ( ^ 

\ M © 

The kinetic energy is given by 



(/ige+0.76/ime). (4) 



E k = l -M W vv 2 = 3.6 x 10 50 erg s" 1 



(5) 



where we have set the v = 6000 km s -1 . For the 
binding energy, we construct a series of WD models 
with central temperature T = 10 7 K using the Models 
for Experiments in Stellar Astrophysics (MESA) code 
(jPaxton et al.1 1201 If ) . We then calculate the binding 
energy as a function of Mwd within the range 0.7 < 
Mwd < 1.4 M© (i.e., the expected range of masses for 
a sub-Chandrasekhar double-detonation model). Mod- 
els of sub-Chandrasekhar-mass WD explosions generally 
find t hat the bulk of IGEs synthesized is in the form of 
56 Ni (|Sim et alJ[20T0l : IWooslev fc Kase^[20TTh . Based 
on our calculated 56 Ni mass in fl4.ll we then have 



/lGE = 0.17 
/iME = 1 



( 



WD 



/iGE- 



and 



(6) 
(7) 



Over the range of WD masses, we find that the nu- 
clear energy released is larger than the combined bind- 
ing energy and kinetic energy. This ensures that there 
is enough energy to unbind the WD and power the ex- 
plosion based on the 56 Ni mass, but it is unclear what 
happens to the excess nuclear energy. Realistically, some 
of the nuclear energy should go into heating the ejecta. 
If some fraction of the WD was left unburned, then E mc 
would decrease. However, our evidence for unburned car- 
bon (see fl3.4p is questionable. Given the uncertainty in 
our assumptions, we cannot place a strong constraint on 
the mass of the progenitor. 

4.3. Late-Time Decay 

We build a simple toy model to describe the luminos- 
ity due to the deposition of energy from the radioactive 
decay of 56 Ni and 56 Co to study the behavior of the 
bolometric light curve. The energy initially deposited 
into the ejecta is from the thermalization of 7-rays emit- 
ted by radioactive decay of 56 Ni — s- 56 Co with an e-folding 
time of 8.8 d. By ~ 20 d, the energy emitted by the de- 
cay of 56 Co— > 5e Fe exceeds that of 56 Ni (e-folding time 
of 111.3 d). The decay of 56 Co can proceed via either 
electron capture emitting a spectrum of 7-rays (96.5% of 
emitted energy) or b eta decay releasin g a positron (3.5% 
of emitted energy) ([Nadvoz hin 1994). The luminosity 



due to the radioactive decay can be written as 

Lrad=£Ni, 7 [l ~e- T ] + (8) 
ico,c+ + Lco,y[l — e _T ], 

where the factor (1— e~ r ) is the fraction of 7-rays trapped 
in the expanding ejecta with an optical depth r. We 
assume that the positrons are fully trapped within the 
ejecta and deposit their kinetic energy instantaneously. 
Note that r is a function of time, decreasing as the ejecta 
expan d h omologously. Following ISollerman et al.l (2002, 
120041) and lStritzinger et al.l (|2006arT 7 we write the optical 
depth to 7-rays in expanding ejecta as 



(9) 



where to is the fiducial time relative to explosion when 
the ejecta become optically thin to 7-rays. 

Solving the first-order differential equations that de- 
scribe the parent-daughter relationship between 56 Ni 
and 56 Co, we write the radioactive luminosity as 



;ad " 



; ANiQNi,7^V"Ni,oe 



>Nii^ 

iV Ni ,o(e- Ac 



(10) 



x [Qco : o+ + Qco, 7 (l ~ e~ r )] , 

where Am and Ac are the inverse e-folding times for 
56 Ni and 56 Co, respectively; A^o is the initial amount 
of 56 Ni synthesized in the explosion; QNi.7(l-75 MeV) 
is the energy yielded by each 56 Ni — > 5e Co decay; and 
<2co, 7 (3.61 MeV) and Q COiC +(0.12 MeV) are the energy 
yielded per 56 Co — > 56 Fe decay via electron capture and 
beta decay, respectively. 

In Figure [121 we plot the bolometric light curve of 
SN 2002es (see gT]for details), as well as the UVOIR 
bolome tric light curves of SN 2005cf f rom iWang et al " 



and SN 1999by using data from lGarnavich et al 
(|2004T ). We also plot models L rad (to = 40 d) and 
L la ,d{to = 00 d) assu ming A/n; = 0.17 M Q (found using 
Arnett's law in fl4.ll) . The model with to — 00 represents 



the case in which the 7-rays are completely trapped and 
the luminosity decays according to the 56 Co decay rate 
(0.01 mag d" 1 ). 

The bolometric light curve of SN 2002es shows a sur- 
prising drop in luminosity at t > +30 d. The decline 
is significantly faster than the cobalt decay rate and the 
decline rate for SN 2005cf. Similar results were found 
for t he late-time decline in individual optical bands (see 
i )3.ip . Even when accounting for a 7-ray optical depth 
that decreases with time, we are unable to find a model 
that reasonably matches the bolometric light-curve de- 
cline rate of SN 2002es. Our best match is to = 40 d, 
which adequately describes SN 2002es within the region 
< t < +30 d but fades too slowly at t > +30 d. 

The decay of the bolometric light curve in the range 
+30 < t < +70 d appears incompatible with the gradual 
escape of 7-rays through homologously expanding ejecta 
and implies that the ejecta become optically thin to 7- 
rays very rapidly. If the ejecta are optically thin to 7- 
rays, then the luminosity should be powered by the ther- 
malized kinetic energy of positrons, assuming some frac- 
tion of the positrons are trapped in the ejecta. We esti- 
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mate the nickel mass required to power the light curve by 
positrons at our last photometry epoch (t = 89.3±3 d af- 
ter explosion) using Equation [TOl and setting r = (i.e., 
the ejecta are optically thin to 7-rays). For complete 
positron trapping, we estimate Mni = 0.05 ± 0.02 M Q is 
required to power the luminosity at this phase; this is a 
factor of 3 less than our previous estimate using Arnett's 
law at maximum light. A corollary to the assumption 
of complete positron trapping is that the light-curve de- 
cay should follow the cobalt decay rate, which is clearly 
not seen. If we allow for only partial trapping, we need 
a positron trapping fraction of ^0.3 to match our Mm 
from Arnett's law. 

The unexpectedly fast decline calls into question 
whether SN 2002es was powered by 56 Ni decay. Here 
we discuss possible ways to reconcile a thermonuclear 
SN with the measured decline rate. 

Dust formation around the SN would lead to a drop 
in optical flux as high-energy photons are reprocessed to 
longer IR wavelengt hs. In models o f dus t formation in 
the ejecta of SNe la, iNozawa et al.l (|2011l) find that the 
conditions required to produce dust occur 100-300 d af- 
ter maximum light. However, due to the low densities 
in SN la ejecta, the dust grains are small (< 0.01 /xm), 
and the expected IR emission associated with dust has 
not yet been detected in late-time observations of SNe la 
(|Gerardv et al.ll2007l ). We also do not find evidence of 
an increase in the red continuum in our spectra or asym- 
metries in line features in our late-time spe ctra as have 
been seen in other SNe that formed dust (jSmith et al.1 
2008). Our last epoch of photometry taken +73 d after 
maximum light gives B — R k, 0.5 mag, while ph otom- 
etry of SN 2005cf taken from I Wang et all (|2009b) gives 
B - R « 0.7 mag. SN 2002es is bluer than SN 2005cf, 
which is opposite the effect we would expect if dust were 
f acilitating the ra pid fading. 

lAxelrodl (|198C0 predicts that SNe la should undergo 
an "infrared catastrophe" (IRC) at late times, once the 
temperature drops below a critical threshold (~ 1000 
K). Models of the IRC predict that a thermal instabil- 
ity shifts the bulk of emission from the optical to fine- 
structure transitions of iron in the IR about 500-700 d 
after maximum light. Without IR data, we are unable 
to determine whether the IR flux increases as the optical 
flux decreases, although the onset of the IRC at such an 
early phase is certainly unexpected. Our last spectrum 
of SN 2002es taken 70 d after maximum light shows per- 
mitted lines in absorption, indicating that the ejecta have 
not yet become nebular and the temperature has proba- 
bly not dropped sufficiently low to facilitate the onset of 
the IRC. The IRC has never been detec ted in late-time 
observ ations of other SNe la, although iLeloudas et al.1 
(2009) evoke the possibility of the IRC occurring locally 
in clumpy ejecta to explain the missing flux in the late- 
time light curve of SN 2003hv. 

4.4. A Pure Explosion Model 

Another possible model for the evolution of the light 
curve of SN 2002es is that the energy deposited into the 
ejecta is derived from the explosion of the progenitor, 
and there i s no subsequent heati ng. Such a model was 
explored by iKasliwal et al.l (|2010l ) to explain the rapidly 
evolving SN 2010X. The explosion energy (Eq) is de- 
posited instantaneously into the ejecta. At maximum 



light, L max ~ Eo/td(0), where id(0) is the initial diffu- 
sion time for photons through the ejecta. The diffusion 
time is td(0) oc M e jK,/Ro, where Rq is the initial radius 
of the SN. Assuming k = 0.1 cm 2 g -1 for Fe-rich ejecta 
(|Pinto fc Eastmad 1200011 ) . L max = 4.0 x 10 42 erg s"\ 
and M e j = 0.66 M e , we estimate an initial progenitor 
radius of R « 10 12 cm. 

This value i s sim ilar to the estimate presented by 
iKasliwafet al. (201(J, in which they argue that such a ra- 
dius would require a progenitor with an extended hydro- 
gen envelope. Based on the absence of hydrogen in spec- 
tra of SN 2010X they reject this hypothesis. Similarly, 
given the lack of hydrogen in our spectra of SN 2002es, 
we also find a pure explosion an unlikely mechanism to 
power the luminosity of SN 2002es. 

4.5. A Core-Collapse SN? 

If SN 2002es was not powered by 56 Ni decay, it is 
reasonable to investigate core-collapse mechanisms that 
could explain the properties of SN 2002es. However, the 
lack of hydrogen emission in optical spectra of SN 2002es 
rules out a massive star with a large hydrogen envelope as 
a progenitor, leaving the possibility that SN 2002es could 
be a SN Ic. The star-formation history of the host galaxy 
strongly favors an older stellar populatio n, making the 
proge nitors of SNe Ibc unlikely as well (jLeaman et al.1 
l20liK If SN 2002es is a core-collapse event, it would 
have to be the result of a low-mass star following an 
a typical evolut i onary path to a SN Ic. 

iPerets et all (|2010l ) discussed data on SN 2005E, a 
low-luminosity SN lb (based on the absence of hydro- 
gen and the presence of helium) with calcium-rich ejecta 
that exploded in the outskirts of its early-type host 
galaxy. Although the observed properties of SN 2005E 
and SN 2002es are different (e.g., SN 2002es is much 
brighter, has slower expansion velocities, and no de- 
tectable helium) , SN 2005E and other objects like it may 
be core-collaps e SNe c onnected to old stellar populations. 
IKasliwal et al.l (|2011[ ) studied the observed properties of 
the emerging class of Ca-rich objects and found that they 
could not be explained by conventional core-collapse or 
thermonuclear explosions. 

Another possible atypical core-collapse object was 
SN 2008ha, a low-luminosity Type I SN (based on 
the absence of hydrogen) that peaked at —14.2 mag 
with ext remely low ejecta velocities (~ 4000-500 
km s- 1 : IValenti et al.l 120091 iFolev et al.l 12001 [2010a]). 
iValenti et al.l <|2009T ) argued that SN 2008ha and the fam- 
ily of SN 2002cx-like objects may be the result of the core 
collapse of hydrogen-poor, low-mass stars. Their analysis 
of SN 2008ha spectra showed little evidence for the IMEs 
commonly associated with the byproducts of a thermonu- 
clear explosion and an absence of forbidden iron lines at 
late times. Coupled with the low luminosity and low 
ejecta velocities, the authors interpret SN 2008ha as a 
core-collapse event that produced little 56 Ni. However, 
I Foley et al.l (|2010aD presented an early-time spectrum of 
SN 2008ha that clearly exhibits IMEs such as silicon, sul- 
fur, and carbon, arguing in favor of the thermonuclear 
explosion of a C/O WD. 

4.6. SN 1999bh: A SN 2002es-Like Object 
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TABLE 5 

Log of Optical Spectral Observations of SN 1999bh 



UT Date 


Phase 3 (d) 


Telescope /Instrument 


Exp. Time (s) 


Observer 1 " 


1999 Apr. 09.3 


+5 


FLWO/FAST 


1200 


MC 


1999 Apr. 24.4 


+20 


Lick/Kast 


1800 


AF, WL 



a Rest-frame days relative to the date of B ma ; 

rounded to the nearest day. 

b AF = A. Filippenko, MC = M. Calkins, WL 



1999 Apr. 03.3 (JD 2,451,271.8), 
W. Li. 
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Fig. 13.— BVRI light curves of SN 1999bh (filled circles) in com- 
parison to SN 2002es (solid lines). The light curves of both objects 
have been shifted relative to the time of -Bmax and peak magni- 
tude. The light curves of SN 1999bh display a striking similarity 
to those of SN 2002es. Both sets of light curves are particularly 
broad, despite being subluminous and spectroscopically similar to 
those of SN 1991bg. 

Combing thro ugh our photometric 

( Ganeshalingam et al.l 120101) and spectral 
( Silverman et al .1 I2012bl) databases, we recognized 
that SN 1999bh shares many of the same properties 
as SN 2002es. Although our photometric and spectro- 
scopic coverage of SN 1999bh is not as extensive as 
our SN 2002es dataset, the available data provide a 
compelling case to link the two objects. 

BVRI photometry of S N 1999bh are taken from 
iGanes halinga m et al.l (|2010f ). Spectra of the object were 
obtained using the Kast dual spectrograph on the Shane 
3-m telescope at Lick Observatory and the FAST spectro- 
graph mounted on FLWO at Mount Hopkins. A journal 
of our observations is available in Table [SJ The spectra 
were reduced using the techniques described in £ 12.21 

SN 1999bh was discovered in NGC 3435 by[Q (fl999h as 



Fig. 14. — Color curves of SN 1999bh (solid circles) compared to 
SN 2002es (solid lines). SN 1999bh has been corrected for E(B - 
V)mw = 0.015 mag and E(B — V\ ost = 0.48 mag. The color 
curves of SN 2002es, corrected for E{B - V)mw = 0.183 mag, 
provide an excellent match for SN 1999bh. 

part of LOSS on 1999 March 29. 2. Subsequent s p ectro - 
scopic follow-up observations by lAldering et al.l (|1999D 
classified the object as a SN la near maximum light on 
1999 Apr. 02. The authors estimate a redshift of 0.028, 
likely based on SN features. The redshift of NGC 3435 
listed in NED is Zhdio = 0.0172. We measure a redshift 
of ^hciio = 0.0168 from narrow Ha + [N II] lines from 
the host galaxy in a spectrum of SN 1999bh taken on 
1999 Apr. 24 using the Kast dual spectrograph mounted 
on the Shane 3-m telescope. We adopt this redshift 
for the remainder of our analysis. The luminosity dis- 
tance to NGC 3435 calculated using zcmb = 0.0172 is 
d L = 70.8 ± 5.3 Mpc. At this redshift, SN 1999bh ex- 
ploded at a projected distance of 3.5 kpc from the nucleus 
of NGC 3435. 

Figure [H shows the BVRI light curves of SN 1999bh; 
we see that those of SN 2002es provide an excellent match 
in all bands. Both objects exhibit broad light curves de- 
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Fig. 15.— Spectra of SN 1999bh (red) in comparison to SN 2002cs 
(black). The spectra have been corrected for the effects of host- 
galaxy recession. Spectra of SN 1999bh have been corrected for 
E(B-V) = 0.01 mag due to Milky Way extinction and E(B-V) = 
0.48 mag from host-galaxy extinction using a CCM reddening law 
with Ry = 3.1. Spectra of SN 2002es have been corrected for 
E(B - V) = 0.183 mag due to Milky Way extinction. 

spite being subluminous, and lack a prominent shoulder 
in the R or I bands. 

Using SN 2002es as a template, we measure the date 
of maximum light to be JD 2,451,271.8 ± 1 d (1999 
Apr. 3.6). After correcting for E(B — V) = 0.015 mag 
from Milky Way ex tinction using the dust maps of 
ISchleeel et all (fl99l . we measure -B ma x = 18.63 ± 0.06 
mag. Due to the faintness of the object, we were un- 
able to follow SN 1999bh sufficiently long to measure the 
light-curve decay rates at t > +30 d. 

Based on the presence of Na I D absorption at the red- 
shift of NGC 3435 in spectra of SN 1999bh, we can be 
certain there is some amount of extinction due to the host 
galaxy. We estimate an equivalent width of 0.8 ± 0.2 A 
of Na I D absorption at the redshift of the host galaxy 
from our spectrum taken on 1999 Apr. 09. Translating 
Na I D EW measurements to an inferred E(B — V) red- 
dening gives 0.1 < E(B — V) < 0.2 mag depending on 
wheth er we use the rela ti onshi p g iven by iBarbon et al.l 
(1990) or iTuratto et alj (120031) . iBlondin et all (|2009f ) 
and IPoznanski et al.l (|201ll ). however, have shown that 
while there is a positive correlation between the presence 
of Na I D absorption and host-galaxy extinction, Na I D 
absorption is not a strong predictor for the amount of 
extinction. 

Instead, we estimate the amount of host-galaxy extinc- 
tion by matching the B — V color of SN 1999bh to that 
of SN 2002es. This assumes that any difference in colors 



is purely associated with host-galaxy extinction and not 
intrinsic differences between the colors of the two objects. 
We find that E(B - V) host = 0.48 ± 0.07 mag. 

We also estimate the host-galaxy reddening by match- 
ing our SN 1999bh spectra to corresponding SN 2002es 
spectra. We perform a fit to match the S EP of the two 
objec ts applying a CCM reddening law (|Cardelli et al.1 
Il989f ) with R v = 3.1 to deredden our SN 1999bh spectra. 
We obtain a best fit of E{B ~ V) host = 0.59 ± 0.01 mag 
(statistical error only) matching our +5 d spectrum of 
SN 1999bh to our +3 d spectrum of SN 2002es, and 
E(B-V)hoat = 0.49±0.01 mag matching our +20 d spec- 
trum of SN 1999bh to our +13 d spectrum of SN 2002es. 
These values are consistent with what we derived above 
using the B — V color at maximum light. We adopt 
E(B - VQhost = 0.48 ± 0.07 mag as the host-galaxy ex- 
tinction. 

In Figure [14] we show the color curves of SN 1999bh 
corrected for Milky Way and host-galaxy reddening. The 
curves of SN 2002es provide an excellent match for the 
color evolution of SN 1999bh. The B-V and V-R colors 
agree almost perfectly. The V — I color of SN 1999bh is 
bluer than expected, compared to that of SN 2002es. 

After correcting for Milky Way and host-galaxy extinc- 
tion, the absolute magnitude at peak of SN 1999bh was 
Mb = —17.71 ± 0.27 mag, almost the same luminosity 
as SN 2002es. We estimate Ami 5 (B) = 1.24 + 0.10 mag. 
Having a light-curve width that is comparable to that of 
SN 2002es, SN 1999bh is another example of a 5cr outlier 
in the Phillips relation as shown in Figure [3] 

In Figure [15] we illustrate our two epochs of spec- 
troscopy compared to SN 2002es at a comparable phase. 
In the top panel, we show our spectra around maximum 
light; SN 1999bh lacks hydrogen and exhibits Si II A6355, 
confirming its classification as a SN la. In addition, there 
is strong Ti II absorption, linking it to the SN 1991bg 
subclass. The ejecta velocity measured from the mini- 
mum blucshift of Si II A6355 is ~6000 km s _1 , similar 
to that of SN 2002es. In the bottom panel, we show 
our spectra from roughly 3 weeks after maximum light. 
Again, we see striking similarities between the two ob- 
jects. We note that the narrow feature at the wavelength 
of Ha in both spectra is from the host galaxy. 

The host of SN 1999bh is classified by 
van den Bergh e~tal\ (2002) as an Sb galaxy. Given the 
narrow Ha, [N II], and Na I D features superimposed on 
our spectra of SN 1999bh, it is likely that SN 1999bh lies 
along the line-of-sight of a star-forming region in NGC 
3435. However, it is not possible to determine whether 
SN 1999bh is actually associated with the star-forming 
region or if SN 1999bh is behind it. Unlike the case 
for SN 2002es, the host galaxy of SN 1999bh does not 
help constrain the stellar population associated with 
SN 1999bh. 

4.7. Rate 

Both SN 2002es and S N 1999bh were included in 
the LO SS SN rate study (|Leaman et all 120111 : iLi et all 
1201 lalf o) in the luminosity function (LF) subsample used 
to calculate the volumetric rate for different SN la sub- 
types. The LF subsample is considered to re present a 
comple te sample of 74 SNe la within 80 Mpc. In ILi et all 
(|2011bf) . both objects are classified as SN 2002cx-like 
SNe. Within a fixed volume, SN 2002es-like objects 
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TABLE 6 

Comparison of Properties of Subluminous SN 1991bg-Like Objects 



SN 


M B (mag) 


Am 15 (B) (mag) 


u p hot a (km s 1 ) 


n(8i) 


Source 


SN 1991bg 


-16.60 ±0.03 


1.93 ±0.10 


10,100 


0.62 ±0.05 


Taubcnbcrgcr ct al. (2008) 


SN 2006bt 


-18.83 ±0.06 


1.09 ±0.06 


10,500 


0.44 ± 0.05 


Foley et al. f2010bT" 


PTF 09dav 


-15.33 ±0.08 


1.87 ±0.06 


6100 


0.35 ± 0.05 


Sullivan et al.(2011b) 


PTF lOops 


-17.66 ±0.06 


1.12 ±0.06 


10,000 


0.58 ± 0.06 


Maguire et al. (2011) 


SN 2002cs 


-17.78 ±0.12 


1.28 ±0.04 


6000 


0.55 ± 0.05 


This work 



As measured by the minimum in the absorption feature attributed to Si II A6355. 



should account for ^2.5% of SNe la. It is worth not- 
ing that the reclassification of SN 2002es and SN 1999bh 
decreases the reported volumetric fraction of SN 2002cx- 
like objects from 5.0% to 2.5% as well. We caution that 
our rate calculation is limited by the small number of 
SN 2002es-like objects in the LF sample. 

4.8. Comparison to Other Peculiar SN 1991bg-Like SNe 

The spectroscopic subclass of SN 1991bg-like SNe 
has been shown to have their own form of Phillips 
relation (i.e., correlation between light-curve decline 
and luminosity), similar to what is found among 



and luminos ity), similar to what is round among nor- 
mal SNe la (iPhillips et all 119991: iGarnavich et all (20041 : 



iTaubenberger et al.l 120081 )7 However, there are notable 
exceptions in the literature of SNe that share superficial 
similarities to SN 1991bg, but have unique properties not 
seen in the broader SN 1991bg subclass. Here we will 
compare and contrast SN 2002es with known peculiar 
S N 1991bg-like SNe 

iFolev et al.l (|2010bD found that SN 2006bt had broad, 
slowly evolving light optical curves (/S.mi§(B) = 1.09 ± 
0.06 mag) typical of a luminous event, but spectra more 
similar to those of a cooler, sublumino us event. A spec- 
trum of the host galaxy from the SDSS (jAbazaiian et al.l 
2009) shows no signs of emission lines, indicating no 
recent star formation. Coupled with the position of 
SN 2006bt in the halo of the galaxy, the authors con- 
cluded that the progenitor of SN 2006bt was likely a 
low-mass star. Assuming no host-galaxy extinction, 
the peak absolute magnitude of SN 2006bt was Mb = 
— 18.94 ± 0.06 mag, slightly less luminous than normal 
SNe la, but not a significant outlier in the Phillips re- 
lation. SN 2006bt had a light-curve width similar to 
that of SN 2002es, but is significantly brighter. Un- 
like SN 2002es, the late-time decay in the light curves 
of SN 2006bt are consistent with 56 Ni decay. SN 2006bt 
had faster expansion velocities at maximum light com- 
pared to SN 2002es. 

PTF has been a prolific discovery engine for unique 
transient events. Two subluminous SNe, in particular, 
may b e related to SN 2002 es: PTF 09dav (iSulliyan et al.l 
I2011b[ ) and PTF lOops (|Maguire et all boill ). Both 
events have maximum-light spectra that resemble those 
of SN 1991bg. PTF 09dav had slow expansion ve- 
locities of ~ 6100 km s- 1 (similar to SN 2002es), 
while PTF lOops had more typical velocities of ~ 
10,000 km s" 1 . Photometrically, both PTF lOops and 
SN 2002es have broad light curves despite being sublu- 
minous. Both objects occupy a similar position off of 
the Phillips relation in Figure |31 PTF 09dav has a more 
typical SN 1991bg-like light-curve shape (Amis(i?) = 
1.87 ± 0.06 mag), but was ~ 1.5 mag fainter than typ- 
ical subluminous SNe la (Mb = —15.5 mag). How- 



ever, SN 2002es is the only object where the optical light 
curves plummet after t > +30 d. 

In Table |6l we summarize the photometric and spectro- 
scopic properties of other peculiar SN 1991bg-like SNe in 
comparison to SN 2002es. While SN 2002es shares some 
characteristics with each of these SNe, no single previ- 
ously published object is exactly like SN 2002es. 

4.9. Models 

Explaining the properties of SN 2002es within the con- 
fines of current models is challenging. Based on the 
star-formation hi story of the host galaxy of SN 2002es 
(jNeill et al.l 120091 ) . we expect the progenitor to be from 
a relatively old stellar population, indicating that it was 
an explosion of either a low-mass single star or a WD. 
In £14.51 we examined the possibility that SN 2002es is 
the result of core collapse of a low-mass single star. In 
this section, we explore possible models to explain the 
observed properties of SN 2002es assuming a WD pro- 
genitor. 

The ".la" model has been proposed to explain rapidly 
evolving, subluminous thermonuclear events. They are 
the result of helium accretion onto a WD from a 
double-degenerate white dwarf A M Canum Venatico- 
rum (AM CVn) bina r y system flBildsten et al.l 120071 : 
iShen fc BildstenT l200l iShen et al.l I2010D leading to a 
thermonuclear explosion of the accreted He envelope. 
These events are expected to have both a luminosity and 
a timescale that is l/10t h that of SN e la. C alculations of 
observable properties bv IShen et al.1 (2010) find that the 
rise time of these objects is < 10 d, with a fast decline 
after maximum. SN 2002es, on the other hand, has a 
rather broad light curve, with a rise that is likely longer 
than 10 d based on a first detection by LOSS of 9 days 
before the time of bolometric maximum light. The ejecta 
velocities are also expected to be ~ 10 4 km s _1 , almost 
a factor of 2 larger than w hat is observed in SN 2002es. 

iWooslev &; Kasenl (|201lD study a broader set of one- 
dimensional simulations with models involving a sub- 
Chandrasekhar-mass carbon-oxygen (CO) WD undergo- 
ing helium accretion, including models that explode just 
the helium envelope and models that explode the entire 
star. The authors survey the parameter space of WD 
masses, accretion rates, and initial WD luminosities to 
produce model spectra and light curves for the resulting 
transient. 

Models in which just the helium shell explode (ei- 
ther through detonation or deflagration) have rise times 
< 10 d and Ami^(B) > 2.0 mag. The fast evolution 
of the B band is caused by the small ejecta mass (de- 
pending on the size of the envelope) and the redistri- 
bution of fl ux from the optical to the NIR by IGEs in 
the ejecta (jKasenl 120061 ). Detonations of "cold" WDs 
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(L = 0.01 L Q ) produce spectra lacking IMEs, which are 
clearly seen in our spectra of SN 2002es. One of the more 
promising models from the set of helium-envelope explo- 
sions involved "hot" WDs (L = 1 L Q ). These explosions 
produce a significant amount of IMEs, resulting in spec- 
tra resembling SN 1991bg with ejecta velocities ^9000 
km s _1 . These models produce a very small amount 
of 56 Ni (~10~ 4 M Q ) and are instead powered by 48 Cr, 
producing a light curve with a 3 d rise time and peak ab- 
solute magnitude of —13. A light curve powered by 48 Cr 
and subsequent decay to 48 V (half life r = 16 d) is an 
appealing explanation for the rapid fading of SN 2002es 
at t > +30 d; however, these models evolve too fast and 
have ejecta velocities that are inconsistent with what we 
observe. 

Deflagrations within the helium envelope led to incom- 
plete burning, producing lower ejecta velocities (^4000 
km s _1 ) and broader light curves caused by the increased 
diffusion time in comparison to the detonation model of 
comparable brightness. The light curves of these models 
are also powered by the decay of 48 Cr, producing faint 
transients (Mb ~ — 14 mag). However, the evolution 
of these light curves is still much too fast compared to 
SN 2002es. 

Outcomes from exploding the entire star generally de- 
pended on the initial luminosity of the WD. "Cold" WDs 
required a larger accreted mass to ignite the helium en- 
velope and drive an explosion of the entire star. As a 
result, the "cold" models contained an outer envelope of 
IGEs (synthesized from helium) which acted as a heating 
source to ionize IMEs and reduce their opacity. Model 
spectra of these events lack Si II, S II, and Ca II, all of 
which are present in SN 2002es. The "hot" models, on 
the other hand, explode with a smaller accreted enve- 
lope, allowing IMEs to appear in the spectra. In fact, 
the synthetic spectra for these objects look remarkably 
like those of normal SNe la and even produce a width- 
luminosity relation, similar to the Phillips relation, but 
with a different slope. However, the minimum expansion 
velocities seen in these models are ~1 1,000 km s _1 , much 
h igher than what is observe d for SN 2002es. 

iPakmor et al.l (|2010L 1201 ID find in simulations that the 
merger of nearly equal mass WDs leads to undcrlumi- 
nous explosions similar to SN 1991bg. Their model fol- 
lows the evolution and subsequent explosion of two WDs 
of equal mass (0.89 M Q ). They find explosions that 
have roughly the same kinetic energy as a normal SN la 
(~10 51 erg s _1 ), but lower velocities due to the larger 
ejecta mass. Less 56 Ni (0.1 M Q ) is synthesized due to 
lower densities in the final merged object, resulting in a 
subluminous event. Synthetic light curves of their mod- 
els are similar to those of SN 1991bg, but broader, with 
Ami5(B) = 1.4-1.7 mag. Synthetic spectra show strong 
titanium, as well as the presence of IMEs such as Si II and 
O I at velocities lower than typically seen in SN 1991bg- 
like objects. All of these match characteristics seen in 
SN 2002es. However, their models do not predict the 
fast drop in flux at t > +30 d, which is a key character- 
istic that makes SN 2002es a unique object. 

5. CONCLUSIONS 

SN 2002es is a peculiar, subluminous SN la with a 
unique combination of observables. At maximum light, 
spectra of SN 2002es are similar to the subluminous 



SN 1991bg, indicating a cool photosphere, but with 
ejecta velocity of ~6000 km s _1 . Such slow velocities are 
more characteristic of SN 2002cx-like objects. While also 
subluminous, SN 2002cx had a maximum- light spectrum 
resembling that of SN 1991T, which is characteristic of 
a hot photosphere. 

Photometrically, SN 2002es has a broad light curve 
(Ami$(B) = 1.28 mag), despite being subluminous 
with a peak absolute magnitude of Mb — —17. 78 mag. 
SN 2 02es is a 5cr outlier in the Phillips relation ((Phillips 
119931 : iPhillips et "all Il999f ) used to calibrate the light- 
curve width vs. luminosity relationship for SNe la. The 
R- and /-band light curves are broad and lack the shoul- 
der typically seen in SNe la. 

Quantitative measure ments of spect ral features such as 
the silicon ratio (^(Si): lNugent et al.l[l99l and the Si II 
A6355 velocity gradient are similar to those of SN 1991bg- 
like objects. However, SN 2002es is an outlier in the 
usual relationships, which show strong correlations be- 
tween the se spectral measurements and light-cu rve pa- 
rameters (|Benetti et al.l l2005: Br anch et alll2006l) . 

From Arnett's law, we estimate a synthesized radioac- 
tive nickel mass of 0.17 M Q required to power the light 
curve. However, the bolometric light curve shows an un- 
expected drop in luminosity at t > +30 d. We are un- 
able to fit the bolometric light curve with a toy model 
of the radiated luminosity that accounts for a decrease 
in the 7-ray trapping function as the ejecta expand ho- 
mologously. If SN 2002es is a thermonuclear event, then 
the ejecta became optically thin to 7-rays in an unex- 
pectedly dramatic fashion. Similarly, we are unable to 
explain the rapid decay by invoking du st formation or 
the infrared catastrophe (|Axelrodlll980fl . Alternatively, 
SN 2002cs could be powered by some other mechanism 
not yet understood. 

SN 2002es exploded in the outskirt s of the early-type 
SO galaxy UGC 2708. SED fitting from |NeilieLaI| WM) 
and emission-line diagnostics indicate that UGC 2708 is 
likely a LINER galaxy with no current star formation. 
This points to a old star, likely a WD, as the progenitor 
to SN 2002es. 

Finding a published model that matches the peculiar 
collection of observables found in SN 2002es is particu- 
larly challenging. We have not found a convincing match 
by comparing SN 2002es to iWooslev fc Kasenl (|201lD 
models of WDs undergoing helium accretion from a com- 
panion star. Models of the merging and subsequen t deto- 
natio n of two equal-mass WDs from lPakmor et all (|2010l . 
1201 If ) are promising, but do not reproduce the drop in 
luminosity one month after maximum light. 

Looking through the e xisting LOSS photometry 
(jGane shalingam ct al. 2010) and the Be rkeley Supernova 
la Pro gram spectral database (BSNIP; iSilverman et alJ 
I2012bf ). we identify SN 1999bh as a probable SN 2002es- 
like event. Spectra and photometry of the object, while 
limited, match many of the observed characteristics of 
SN 2002es. Both objects are inc l uded in the LO SS SN 
rate studies (jLeaman eFaI]l20Tll [Li et alJl2^TTallbh . al- 
lowing us to estimate that roughly 3% of SNe la should 
be SN 2002es-like SNe within a fixed volume. 

Ongoing surveys for transient objects, such as PTF 
and Pan-STARRS, are likely to find more objects like 
SN 2002es. PTF has already published results for two 
peculiar subluminous objects which share some similari- 
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ties to SN 2002es, but also have distinct differences. We 
expect that future datasets will shed light on understand- 
ing the bizarre nature of these subluminous objects. 
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